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We review recent work by the ALPHA and UKQCD Collaborations where masses and matrix elements were 
computed in lattice QCD using Schrodinger functional boundary conditions and where the strange quark mass 
was determined in the quenched approximation. We emphasize the general concepts and our strategy for the 
computation of quark masses. 



1. Introduction 

Based on the recent progress in lattice QCD 
concerning non-perturbative renormalization of 
local composite operators and removal of lead- 
ing discretization errors [0^, there is a natural 
interest to calculate hadron masses and matrix 
elements in the continuum limit. Especially, the 
problem of a determination of the light quark 
masses can be addressed with confidence, 
since the complete renormalization is known non- 
perturbatively 

As an alternative to more standard methods, 
we will show that the Schrodinger functional (SF) 
framework allows for a reliable computation of 
spectral quantities with high accuracy in lattice 
QCD. As a necessary prerequisite, we will demon- 
strate explicitly that SF correlation functions are 
dominated by hadron intermediate states at Eu- 
clidean time separations of around 3 fm. 

This technique is then applied to physical ob- 
servables in the meson sector of QCD, in particu- 
lar to light quark masses. We present an overview 
of the work discussed in refs. 0J^- AH results re- 
fer to the fully 0(a) improved theory, but due 
to the rather introductory character all issues re- 
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lated to 0(a) improvement are ignored here. 

Our quenched data stem from four lattices (SF 
boundary conditions) of spatial extent L « 1.5 fm 
and time extension T — 2L with lattice spacings, 
a, from 0.1 fm (/3 = 6.0) to 0.05 fm {/3 = 6.45), 
and with four values of the quark mass each, so 
that all results can be extrapolated first in the 
pseudoscalar mass and finally to the continuum 
limit. For our choice of parameters, chiral per- 
turbation theory predicts negligible finite size ef- 
fects and an explicit investigation Q , which is not 
discussed here, confirms this. 

2. SF correlators at large time separations 

The SF is defined as the QCD partition func- 
tion in a L'^ X T cylinder with periodic boundary 
conditions in three of the four Euclidean dimen- 
sions and Dirichlet boundary conditions in time 
at the hypersurfaces xq = and xq = T [d\. 
Whereas in earlier apphcations, e.g. in refs. |^,3|, 
correlation functions in the SF were mainly con- 
sidered in the perturbative regime (i.e. small ex- 
tensions of the space-time volume), the emphasis 
in this section is on their properties in intermedi- 
ate to large volumes with extensions significantly 
larger than the typical QCD scales of order 1 fm. 

Starting from the quantum mechanical inter- 
pretation of the SF pjlCl], one can derive explicit 
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expressions for the representation of its correla- 
tion functions in terms of intermediate physical 
states 0. Let us discuss an example which is 
directly relevant in section 3. We consider the 
correlation functions 



fx{xo) = -^{X{x)0), 
h = -UO'O), 



(1) 
(2) 



where O is a pseudoscalar field, constructed from 
a p = quark field and a p = antiquark field 
at a:o = : O = Ey,z 0(y)75Cj (z) , be- 
ing flavor indices. Analogously, O' is located at 
the boundary xq = T. One may think of the 
fields Ci ■ ■ • a-s quark fields at the boundary; their 
precise definition is given in In the fol- 

lowing, we shall be interested in /a, defined by 
X{x) = Ai){x) = ipj{x)joj5ipi{x), and /p, given 
hy X{x) = P{x) = ^j{x)j5M^). 

For large separations xq and T — xq the correla- 
tion functions are dominated by the lowest lying 
intermediate states with the appropriate quan- 
tum numbers. These are the pseudoscalar ground 
state and the vacuum, the latter contributing be- 
tween Xq and T in both /a and /p. In the trans- 
fer matrix formalism one may then arrive ||^ at 
the asymptotic behaviour including the first non- 
leading corrections: 



fxixo) 



h 
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(3) 



In leading order an unknown matrix element 
p and the desired matrix elements of the 
Schrodinger picture operators X , associated with 
the fields X , arise. The mass of the 0++ glueball, 
mc, and the gap in the pseudoscalar channel. A, 
enter the corrections with coefficients rfx ? 
pressed through certain matrix elements. Apart 
from the pseudoscalar mass, mps, the formulae 
(^) enable to calculate the decay constant, .Fps, 
and the pseudoscalar coupling, Gps, via 

(0,0|Ao|0,PS) = Fpsmps(2mpsL3)-l/^ (4) 
(0,0|P|0,PS) = Gps(2mpsL3)-V2. (5) 

Note in particular that /i may be used to elimi- 
nate the unknown factor p, the only place where 
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Figure 1. The ratio /a//p for a « 0.07 fni (open 
circles) and a « 0.09 fni (filled circles). 

the divergent renormalizations of the boundary 
quark fields are hidden. A generalization e.g. to 
the vector meson channel is straightforward. 

We turn to the numerical tests of the practica- 
bility of our method. As a representative exam- 
ple, which moreover will play a crucial role in the 
next section, we discuss the ratio 



{1 + (e-="o^, Q-iT-xo)mG^ I 



/a _ "ips-Fps 
/p Gps 

It shows approximate plateaux between « 1 fm 
and T - w Ifm (Fig. p. Their height deter- 
mines TOpsJps/Gps. Before, however, just read- 
ing these off, the magnitude of the "contamina- 
tions" by excited state contributions should be 
assessed. Given the fact that estimates for A 
and mc are available this may be achieved 
by plotting the ratio against the leading correc- 
tions: the slopes in Fig. |2| allow to deduce the 
range tmin < a;o < imax where the excited state 
corrections are below a certain systematic error 
margin, which we chose in turn to be negligible 
compared to the final statistical errors. The ratio 
was then averaged in the range tmin < a^o < tmax- 
This procedure is easily extended to the calcu- 
lation of TOps , Fy>s and Gps separately. 

The numerical efficiency of the computations 
with SF correlation functions is illustrated in 
Fig. H, where we compare our results for the bare 
(improved) pseudoscalar decay constant with 
those found in the literature. To judge the rel- 
ative size of the statistical errors one should note 
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Figure 2. Influence of excited state contributions 
on /a//p at the same parameters as in Fig. 0. 
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Figure 3. Comparison of some of our results 
(ALPHA) with those of the QCDSF ^ and 
APE CoUaborations. Symbols are slightly 
displaced. 

that our statistics is not much higher than the 
statistics present in the other simulations. 

Clearly, the SF applied for extracting hadron 
physics is similar to the method of wall sources, 
but one should note that gauge invariance is kept 
at all stages of the present formulation. Further- 
more, dimensionless non-local fields are used to 
create the boundary states. Dimensional analy- 
sis then tells us that the pre-asymptotic decay of 
SF correlation functions is slow, leaving large and 
precise signals at hadronic length scales of 1-2 fm. 
Since the correlation functions are renormalizable 
by simple factors, this property is indeed indepen- 
dent of the lattice spacing, once one is sufficiently 



close to the continuum limit. 

3. Quark masses: the strategy 

Before we come to the numerical applications of 
the above methods, we want to explain in detail 
our strategy for the computation of light quark 
masses. Their ratios can be computed from chi- 
ral perturbation theory [Q. The currently best 
results read 

Afu/A/d = 0.55 ± 0.04 , Ms/M = 24.4 ± 1.5 (6) 

with M = + Md) |15|. 

Despite this success, there is substantial work 
to be done using lattice QCD. 

• The applicability of chiral perturbation the- 
ory needs to be checked. This concerns the 
practical question in how far the lowest or- 
ders dominate the full result. 

• The parameters in the chiral Lagrangian (at 
a given order in the expansion) can not be 
inferred with great precision from experi- 
mental data alone. Their determination can 
be improved significantly by help of lattice 
QCD results. 

• In particular, there is one parameter in the 
chiral Lagrangian which is impossible to de- 
termine from experimental data. This is the 
overall scale of the quark masses, which is 
only defined once the connection with the 
fundamental theory, QCD, is made. 

An important point to realize is that all of the 
above problems can be dealt with by working with 
unphysical quark masses (of course they should 
not be too large). For the first two problems it 
is in fact essential to explore a range of quark 
masses, and also for the last problem, which we 
address here, this is of significant advantage. 

The above considerations lead us to determine 
an implicitly defined reference quark mass M^of, 



m|g(Mi.ef)ro = 1.5736 = (mxro)^ , 



(7) 



where mpg(Af) is the pseudoscalar meson mass 
as a function of the quark mass (mass-degenerate 
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quarks). Chiral perturbation theory in full QCD 
relates M^ef to the other Hght quark masses viz. 



(8) 



Further evidence for this relation is given in ^J, 
using the results of quenched lattice QCD. 

What was said above is valid literally in mass- 
independent renormalization schemes, where 
quark masses are renormalized with a flavor in- 
dependent factor. The resulting running quark 
masses, m{^J.), are however scheme- and scale- 
dependent quantities. It is advantageous to com- 
pute directly the renormalization group invariant 
(RGI) quark masses which are pure numbers 
and do not depend on the scheme. In terms of 
fn{fi) they are given as 



M 



lim 



{{2b„g\^l)r''"/''>"fn{p)] , (9) 



bo = 11/(4^)', da = 8/(4^)2 



(10) 



For 0(a) improved quenched QCD, the ALPHA 
Collaboration has determined the renormaliza- 
tion factor Zm relating the bare current quark 
masses m, defined by the PCAC relation, to the 
RGI masses [^,||: 



M = Zy\m, . 



(11) 



Applying the PCAC relation to the vacuum- 
to-pseudoscalar matrix elements then results di- 
rectly in our central relation. 



2roM„f = Zm ^'""'"^""' '^"'' X 1.5736. 
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R 



4. Results 



Fps/Gps ■ 



(12) 
(13) 



With the methods of section 2, the ratio R/a 
and the meson mass nipsa can be computed ac- 
curately as a function of the bare quark mass and 
the bare coupling. Also the scale ro/a is known 
|l6|| . A mild extrapolation, shown in Fig. 0, then 



yields R/a at the point 



rgmpg 



1.5736.P Next 



the full combination roMj-ef (see eqs. (|7|,|8|)) is 

^ The reason for not computing directly at or below the 
quark mass corresponding to r^mp^ = 1.5736 was that 
this is the region where unphysical zero modes of the 0(a) 
improved Dirac operator are relevant ("exceptional config- 
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Figure 4. Mass dependence and extrapolations 
at the two smallest values of the lattice spacing. 

extrapolated to the continuum limit in Fig. |^a. 
Since this latter extrapolation involves a signifi- 
cant slope we discarded the point furthest away 
from the continuum. This is a safeguard against 
higher order lattice spacing effects. 

Furthermore, the entire analysis was repeated 
for Mref in units of the Kaon decay constant 
In this case the lattice spacing dependence 
is weaker (Fig. ||b). The two results. 



2 roMref = 0.36(1) 

rp — 0.5 fm 

2AW(Fk)r = 0.87(3) 

(Fk)r = 160McV 



2M„t = 143(5) MeV, 



2Mrcf = 140(5) MeV, 



agree better than one would expect them to do in 
the quenched approximation. 

It is premature to conclude that this approx- 
imation yields a nearly unique result. Rather, 
the assignment of physical units is ambiguous. 
One may estimate [|j by help of the recent results 
of the CP-PACS Collaboration @ that roughly 
10% larger numbers would be obtained, if the 
scale ro were replaced by one of the masses of 
the stable light hadrons. This represents a typi- 
cal ambiguity of the quenched approximation. 



urations"). It is now understood that this technical prob- 
lem can be resolved by introducing tha^quark mass into 
the lattice theory in a particular way pj] . However, even 
without entering this region, our quark mass points are 
close enough to perform a safe extrapolation. 
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Figure 5. Continuum extrapolations of A/rcf in 
units of To (a) and (Fk)r (b). Full symbols are 
the a — limits. The (dashed) fit functions are 
continued beyond their range towards larger a. 

Other interesting results after extrapolation to 
the continuum limit are 



{?'o(-Pps)r}^ 



= 0.415(9) 



compared to 0.5 fm x (^k)R = 0.405(5) with the 
experimental value of {Fk)r, and 



{^OTOvcctorlj 



= 1.5736 



2.39(7) 



to be confronted with 0.5 fm x tok* = 2.26. In 
the latter comparison one should remember that 
the K* meson is unstable even in absence of elec- 
troweak interactions. Its width, Fk*, amounts to 
0.5 fm X Fk* = 0.13 . 

5. Discussion 

In addition to the - by now well known - ap- 
plications to renormalization problems in QCD, 
the Schrodinger functional has been shown to be 
useful for the computation of low energy matrix 
elements. For the particular cases of vacuum-to- 
pseudoscalar matrix elements, our results are sig- 
nificantly more precise than those from standard 
methods. This implies a precise computation 
(in the continuum limit) of the renormalization 



group invariant quark mass Afref defined above. 
MS masses for finite renormalization scales ^ are 
obtained through perturbative conversion factors 
known up to 4-loop precision. A typical result is 



ms(2GeV) = 97(4) MeV , 



where the uncertainty in A^^ = 238(19) MeV g, 
entering the relation of the running quark masses 
in the MS scheme to the RGI masses, has been 
accounted for and the quark mass ratios from full 
QCD chiral perturbation theory were used. 

One should keep in mind that there is an in- 
trinsic ambiguity when assigning physical units 
in the quenched approximation. We estimated 
this to be of order 10%. 
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